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' We study physical properties of conformal initial value data for single and binary black hole 

' configurations obtained using conformal-imaging and conformal-puncture methods. We investigate 

how the total mass Mtot of a dataset with two black holes depends on the configuration of linear or 
' angular momentum and separation of the holes. The asymptotic behavior of Mtot with increasing 

separation allows us to make conclusions about an unphysical "junk" gravitation field introduced 
in the solutions by the conformal approaches. We also calculate the spatial distribution of scalar 
invariants of the Riemann tensor which determine the gravitational tidal forces. For single black 
hole configurations, these are compared to known analytical solutions. Spatial distribution of the 
invariants allows us to make certain conclusions about the local distribution of the additional field 
in the numerical datasets. 
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I. INTRODUCTION 

The problem of initial conditions (initial Cauchy data) for the integration of the evolution of colliding black holes 
is an important problem of the numerical general relativity and has attracted attention of many researchers. For a 



I . recent review see H. One of the approaches to the problem is a conformal-imaging method proposed in ||^, g], |J] 

^jpj' and developed in [Q, Another is the conformal-puncture method |^]. It is known that initial black hole data 

• • , constructed using both these methods contains some non- vanishing dynamical component |jl|], an unphysical, junk 

. ^ ' gravitational field. Without the integration of the initial data in time it is impossible to make exact conclusions about 

k>( ] the character and the amount of this unphysical field. Still it is possible to make some conclusions about the junk 

; I ' field using only the initial data. For example, one can numerically calculate global characteristics of the initial data 

' such as total energy, and then compare it with known analytical or approximate solutions |^ , Q , . 

In our previous paper we developed adaptive mesh refinement approach to the construction of initial data for black hole 
collisions on high resolution Cartesian meshes ||^. The method allows us to compute initial data with high accuracy 
both near and very far away the from black holes. The goal of this paper is to use this method to systematically 
analyze the physical properties of black hole initial data for a wide range of colliding black hole configurations. To 
analyze the junk fields present in the data, we will use two different approaches. First, we calculate and compare 
global characteristics of the configurations. In addition, we use a new approach which consists of calculating and 
comparing four local scalar invariants of the gravitational field ( section 92) . These scalar invariants completely 
determine local properties (tidal forces) of a gravitational field. 

In the next section |^ we briefly describe our adaptive mesh refinement method of calculating initial black hole data. 



We start our discussion of local and global characterization of the solutions in section III . Results for single black 



hole configurations with angular or linear momenta are presented in section [V . For these black holes we can compare 
numerical initial date with exact analytical solutions. In section ^ we present results for two-black hole configurations 
with different orientations of linear and angular momenta and with different separations between the black holes. 
All calculations are performed using both conformal-imaging and conformal-puncture methods, and the results are 



compared. Conclusions are presented in section VI 



1 



II. CONSTRAINT EQUATIONS AND THE METHOD OF SOLUTION 



The ADM or 3+1 formulation of the equations of general relativity works with the metric gf^- and extrinsic curvature 

Kf^ of three-dimensional spacelike hypersurfaces embedded in the four-dimensional space-time, where i,j = 1,2,3, 

and the superscript ph denotes the physical space. On the initial hypersurface, gf^ and Kfj^ must satisfy the constraint 
equations The conformal approach assumes that the metric is conformally flat, 

Sf/ = (t>^9^j , (1) 

where gij is the metric of a background flat space. This conformal transformation induces the corresponding trans- 
formation of the extrinsic curvature 

Kfl^ = r'K,, . (2) 

With the additional assumption of 

trK = , (3) 



the energy and momentum constraints are 



+ Icjy-^'K.jK'' = 0, (4) 



and 



0, 



(5) 



respectively, where and Dj are the Laplacian and covariant derivative in flat space. 

A solution to and (|^) for two black holes can be specified by six parameters (three parameters per black hole). 
These are the mass parameter Ms, linear momentum parameter Pg and angular momentum parameters Ss, where 
6 = 1,2 is the index of a black hole. In terms of these parameters, the solution to the momentum constraint for two 
black holes is 

K,,{r)=Kl^(r)+K^^{r) , (6) 



where 



3 1 

^'r(i") = 9 H 152 {Ps.^ns,J + P5.Jns,^ - {g^J - ns,^ns,j) Ps,knl) (7) 




In and (0) the comma separates the index of black hole from the coordinate component indices, Rs = Ms/2 is the 
black hole throat radius, and = {r — rs)/\r — rs\ is the unit vector directed from the center of the 5-th. black hole 
to the point r. We work in units where G = 1, c = 1 which means that Ps is in mass units and Ss is in mass^ 
units. 

In the conformal-imaging method, we obtain an inversion-symmetric solution to by applying an infinite series of 
mirror operators to (^), as described in [^. Note, that before applying the mirror operators to K^'J^^ , this term must 
be divided by 2 since the image operators will double its value. The series converges rapidly, and in practice only a few 
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terms are taken. After the isometric solution for Kij is found, (Q) must be solved subject to the isometry boundary 
condition at the black hole throats 

n^A^^-^. (9) 

In the conformal-puncture method, expression is used without mirror imaging, and the energy constraint is solved 
on with a puncture at the center of each hole. Brandt and Britgmann |0] proved that it suffices to solve the energy 
constraint everywhere on without any points cut out. Thus, this method avoids inner boundaries. 

The outer boundary condition in both cases is </) ^ 1 at infinity. This boundary condition is represented by ||] 

(10) 

ar r 

where r is the distance from the center of the computational domain to the boundary. 

As described in g, (|) and (|) are solved on an adaptive Cartesian mesh. The mesh is refined on the level of individual 
cells using a fully threaded tree structure Computational cells can be characterized by the level of cells I in the 
tree. A cell size A; is related to the cell level as A; = where L is the size of the computational domain. In 

this paper, L = 1024 and the minimum and maximum levels of refinement is Imin = 5 and Imax = 16, respectively. 
That is, the coarsest cells had a size Amm = 32 and the finest cell size was Amax = 1/64. The unit of length is the 
throat radius Rs = Ms/2 of the holes (in this paper, all holes have the same mass Ms = 2 and thus the same throat 
radius). The equivalent uniform grid resolution thus corresponded to that of a 65536^ uniform grid. Mesh was refined 
at lmi7i < I < Imax according to the following refinement criterium (see ||^ for details) 

(11) 




Mesh was refined at levels ^ = 5, 6 if ?7 < e = 0.03 and at Z > 6 if 77 < e = 0.07. In our previous paper it 
was demonstrated that our mesh refinement approach gives a quadratic convergence of solutions with increasing 
numerical resolution. In this paper we additionally checked the accuracy of numerical solutions by varying e, Imax 
and Imin- The accuracy of the results is discussed further in sections |^ and 



III. GLOBAL AND LOCAL CHARACTERISTICS OF INITIAL DATA 

Global physical properties of a system are its total mass Mtot, linear, Ptot, and angular, Stot, momentum For 
conformally flat solutions 

Mtot^~^{ ipT)da\ (12) 



27r J^^^ dx^ 

P^,tot = -^f K.jda^ , (13) 

5,.tot = ^SMj I {x'K^" - x^K"')dak , (14) 

where integration is performed over the surface at infinity. For our case of two black holes, the total linear and angular 
momenta are equal to 

Ptot=Pl+P2, (15) 
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and 

Stot = Si + S2 + (Ci - O) X Pi + (C2 - O) X P2 , (16) 

where Ci, C2 are centers of black hole throats or punctures and O is the center of mass of the system (see [|l^). Note, 
that Pi, P2, Si, S2 do not have direct physical meaning but Ptot and Stot do. 

In addition to global characteristics of a system, we can also consider local invariants. In general case, there are four 
such invariants ( section 92) 

h = R'^'^'^'Rcp^s (17) 
h = R'"''''Rc(3^s (18) 
/3 = Ro.p-,sR^'^''Rf, (19) 
h = Rcp^sR^'^" Rf, (20) 

where Rap^s = ^^aptiuRjs ^^'^ ^aPfii/ is a fully antisymmetric unit tensor in curved coordinates. At every point, in a 
specially selected tetrad, the Riemann tensor can be expressed in terms of these invariants (except of special degenerate 

cases P). Since Riemann tensor uniquely determines tidal forces, these invariants fully characterize physical properties 

1/2 

of the field at every point. The tidal force is oc /i by the order of magnitude. For a Schwarzschild black hole, for 
example, the tidal force is exactly Mjc^r^ = ^JTi 

In order to calculate the Riemann tensor, we need to know the space and time derivatives of the four-dimensional 
metric at each point on the initial slice. This data can be obtained from the time-evolution part of the Einstein 
equations. Since the scalars do not depend on the choice of a coordinate system, we can select zero shift and unit 
lapse to simplify calculations, and to write the derivatives as 



dt 



2K,j , (21) 



^'^9^] ^ _2^^ (22) 
^ = -2( + tvKK.j - 2K,kK)) . (23) 

Here both gij, Kij are the physical metric and extrinsic curvature, and not their counterparts in flat space, and ^^"^ Rij 
is the Ricci tensor of 3-dimensional space which depends only on the spatial metric and it's spatial derivatives. 



IV. SINGLE BLACK HOLES 



A. Rotating single black hole 

We consider a single black hole with a non-zero angular momentum S 7^ and zero linear momentum P = 0, and 
compare numerical conformal solutions with the known analytical solution for a Kerr black hole. The Kerr metric is 

, ■? ( 2Mtotr\ , AMtotra siii^ 6 E , , ^sin^6' , , , 

ds^ = -[1 ^ ) dt^ dtd(f> + —dr^ + Ed9^ + — — d(l>^ , (24) 



where 



E = r^a^ cos^ 9 , (25) 
A = r^- 2Mtotra^ , (26) 
A = {r^ +a^f -a^l\&in^e . (27) 
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The analytical solution (19), (20) is described by two parameters, the total mass of the hole Mtot measured at infinity, 
and the specific angular momentum a = S/Mtot such that a < Mtot- The irreducible mass of a Kerr black hole, i.e., 
the mass of the event horizon of the hole is given by 



\ 



( s 



V Mt, 



(28) 



We have calculated datasets with single rotating black holes using both the imaging and puncture method. In each 
case, we vary the S parameter, but hold the mass parameter of the black hole, M — 2, constant. We have checked 
the accuracy of the solutions by varying Imax-, Imin and e in (pi]). For each dataset we computed Mtot and Stat by 



evaluating (12) and (14) over at sphere with radius r < L/2 ~ 500. We found that the values of Mtot, 3*°* calculated 
at this radius were third digit accurate. The integration in ( p^ and (^4|) is supposed to be over at surface at infinity, 
and we determined Mtot{oo) and Stot{oo) from the corresponding values at finite radii by extrapolation in 1/r to 
1/r = 0. 

Figure |l| shows the extrapolation procedure for a rotating black hole with S — 40 calculated with the imaging method. 
Mtot and Stot have been calculated by integrating over spheres with increasing radii, r. We find a least square fit of 
the data to a second order polynomial and evaluate that polynomial at r = cx3 to obtain Stotioo) and Mtot{oo). For 
rotating black hole, the values of Mtot are practically not changing with r, whereas the values of Stotioo) are slightly 
larger than the values at finite r. For a case S — 40, the extrapolated value Stot{oo) ~ 40 ~ 5, as it should be. The 
asymptotic value of Mtot is in a good agreement with that obtained in ]l3[ |. Figure |^ shows Mtot{oo) as a function of 
S. 

In the conformal-imaging approach, the location of an apparent horizon for a single rotating black hole coincides 
with the throat [|l3|. Thus, in this case it is possible to investigate the properties of the apparent horizon. We have 
computed the area of the horizon Aah, the mass, Mah = \/ ^ah/167r and the polar and equatorial circumference, Cp 
and Ce- The values of Mah agree within three digits with the corresponding values in for S = 10. Figure || shows 
calculated Mtot and Mah as a function of S for a rotating black hole with M = 2. Also shown is the irreducible mass 
Mirr calculated for a Kerr black hole with the corresponding S and Mtot- Figure ^ shows the calculated ratio of the 
polar and equatorial circumferences Cp/Ce as a function of S and compares it with the corresponding ratio for a Kerr 
black hole 0. 

Figure ^ shows that for the same total mass of the configuration Mtot and the same angular momentum S, the mass of 
the apparent horizon in the numerical data, Alah, is less than the irreducible mass of the black hole. This is expected 
since in the conformal approach there is an additional energy/mass related to a junk field outside of the black hole. 
Figure ^ shows that Cp/Ce is close to unity for conformal imaging solutions whereas for Kerr black hole it is less 
than one. This means that the shape of the apparent horizon is substantially closer to spherical symmetry in the 
conformal-imaging solutions than for a Kerr black hole. 

We have also compared the scalar invariants (17) of our numerical solutions with the corresponding analytical values. 
The difference between the values of these invariants of the numerical and and analytical solutions characterize the 
junk gravitational field. We have chosen to compare the value of the invariant in the points that are the same physical 
distance from the apparent horizon along the equatorial axis of the black hole, because this gives us a correspondence 
between the Boyer-Lindquist radial coordinate r of the Kerr metric and the isotropic radial coordinate x of the 
numerical solution. In the numerical solution the apparent horizon coincides with the throat, and in the Kerr 
solution the apparent horizon coincides with the event horizon. To find the Boyer-Lindquist radial coordinate r that 
corresponds to a given isotropic coordinate x we find the physical distance: 

ARphix) = f ^^dx , (29) 
Jr 

where x is the distance to the center of the black hole throat in the background flat space and R is the throat radius. 
The integration is performed numerically using cubic spline quadrature. We normalize this distance: 

Alphix) = ARph-^ (30) 

-IVltot 

and we write down the following implicit equation: 
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Alph{x)Mtot = 



(31) 



and solve it with respect to r using Van Wijngaarden-Dekker-Brent's method. Mtot is the total mass of the numerical 
black hole with mass parameter M = 2R, as well as the total mass of the Kerr black hole that we are comparing to. 
Finally we compute the value of the invariants for the Kerr black hole (which are known analytically) at r{x). 

Figure (H) shows that increasing value of S leads to the decrease of Ii at the black hole horizon for Kerr black holes. 
This is a well understood effect: from we know that the total mass of a black hole and its apparent horizon grow 
with S and we know that in the Kerr solution tidal forces at the event horizon of a Kerr black hole decrease with the 
increase of the mass. We can see that the same tendency holds for conformal black holes (Figure ^ a,b). Figure ^ c 
shows the difference in Ii for the conformal and Kerr solutions. At the boundary of the black hole the absolute value 
of the difference initially decreases with increasing S. However, for large S it begins to decrease due to the decrease 
of tidal forces. Figure ||d shows the relative difference in Ii . The absolute value of the relative difference at the black 
hole boundary increases with increasing S monotonically. Far away from the black hole, the absolute value of the 
difference is also monotonically increases with S, and for large values of S it is significantly larger than that at the 
black hole boundary. Figure || thus shows that there exists a junk gravitational field in the conformal solutions, and 
that the relative amount of this junk field increases outside of the black hole with increasing S. Figure ^ shows the 
relative difference of Ii at large distances. The relative difference is remains large even at very large r ~ 500. 



B. Moving single black hole 



Now consider a single black hole with a non-zero linear momentum P 7^ and zero angular momentum S = 0. Similar 
to rotating black holes, we compute numerical conformal solutions with M — 2 and various P, calculate Mtot and 
Ptot for these solutions by integrating (12) and (13) over spherical surfaces of different radii r, and then extrapolate 
these values to infinity. Figure illustrates the extrapolation procedure for the case M = 2 and P — 40. In the case 
of a boosted black hole, both Mtot and Ptot vary noticeably with r. The extrapolated value of Ptot (00) ~ 40 ~ P, as 
it should be. We compare our numerical solutions with known analytical solutions for a boosted Schwarzschild black 
hole with the same Mjot (00) and Ptot (00). 

The spatial metric of a boosted Schwarzschild black hole can be found by performing a Lorcnz transformation 

t = -f{t + vx) , (32) 
X — ^{x + vt) (33) 

(34) 



of a Schwarzschild metric, where 7 = (1 — u^)^/^. This gives the following spacetime metric 

\/ T ( t V ,T,d ) \ 1 ,T,4 / 



\I/4 



with 



^ = 1 



M\ 
2r) 



M 
2R 



M 
2R 



R= 



(36) 



The analytical solution ( |35| ) is described by the two parameters, the mass of the Schwarzschild black hole M and the 
velocity of the hole v as measured by an observer at infinity. The total mass and velocity of a boosted Schwarzschild 
black hole are 



P2 P 



Mtot = V^^^, ^-Xlj^2-J^^- (37) 
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At the horizon, the metric is exactly the same as for a Schwarzschild black hole (this can be seen from the fact that 
a = at the horizon). Thus, the horizon properties are the same as for a Schwarzschild black hole. Also, since we 
have only performed a coordinate transformation, /i for a boosted Schwarzschild black hole is the same as Ii for a 
Schwarzschild black hole. 

In the conformal-imaging approach, the location of an apparent horizon for a moving black hole coincides with the 
throat [T^ ] and for a Schwarschild black hole the apparent horizon coincides with the event horizon. We have computed 
the area of the horizon Aah, the mass, Mah — yAi/i/lBTr and the polar and equatorial circumference, Cp and Cg. 
Figure || shows calculated Mtot and Mat as a function of P for a moving black hole with M = 2. Also shown is the 
analytical horizon mass Mboosted — \/Mfot-~^ ■ Figure || shows that for the same total mass of the configuration 
Mtot and the same linear momentum P, the mass of the calculated apparent horizon Mah is less than the analytical 
one. This is expected since in the conformal approach there is an additional energy related to a junk field outside 
of the black hole. The polar and equatorial circumferences for a moving black hole are practically the same, i.e. the 
symmetry of the horizon is close to spherical symmetry. 

We have compared the scalar invariants (17) of our numerical solutions with the corresponding analytical values. We 
have chosen to compare the value of the invariant in the points that are the same physical distance from the apparent 
horizon along the axis of motion (the x-axis). This physical distance can be expressed as a function of the isotropic 
coordinate distance x from the black hole center along the axis of motion to the point where we make the comparison. 
Figure ^ compares Ii as a function of x for various values of P. We see that the increase of the total mass with 
increasing P leads to a decrease in the value of Ii at the black hole horizon (Figure ^ a,b). Figure |^ shows the 
difference in Ii for numerical and analytical solutions. At the black hole boundary, the absolute value of the difference 
initially decreases with increasing P. However for large P it begins to decrease due to the decrease of tidal forces. 
Figure |^d shows the relative difference in /i . The absolute value of the relative difference at the black hole boundary 
first increases with increasing P but than starts do decrease. Far away from the black hole, values of Ii for boosted 
black holes are much larger than Ii for conformal black hole solutions (Figure |l^) . Similar to the rotating black hole 
case, Figure |9| show that conformal solutions contain a junk gravitational field. The amount of the junk field increases 
outside of the black hole with increasing P. 



V. BINARY BLACK HOLES 



There is an infinite number of possibilities for choosing the parameters that characterize the initial data for two black 
holes. The problem can be parameterized by 

(r2,ri,d=^— ^,P^5i,P2',5^) (38) 
2ri 

where = and Mg is the mass parameter of the (5'th black hole. Ci and C2 are the coordinates of the center of 

the throats and punctures respectively. We chose to call the distance d = ]^^^ ^ the separation of the hole. Since 
the holes are placed at equal distance from the center of the computational domain, d is the coordinate distance from 
the center of each hole to the center of the domain. The rest of the parameters are the linear/angular momentum 
parameters. 

We restrict ourselves to some specific case. We consider configurations with equal mass black holes. Mi — M2 — 
M = 2, i.e. Ri = R2 = 1. The line connecting the centers of the holes coincides with the x-axis. We only investigate 
configurations where the holes have either linear momenta or angular momenta, but not both. The configurations 
considered are shown schematically on Figure Arrows on this figure show the direction of either linear or angular 
momenta. In all cases when momenta of both black holes are not zero, they have equal absolute values except when 
explicitly stated otherwise. 

First we mention the following general property of our initial data. Let the linear momentum parameters of the 
problem are Pi and P2, and the angular momentum parameters for both holes are zero. Then the metric that is a 
solution to the constraint equations, ^ and (||) with these parameters is the same as the metric of that is the solution 
to the constraint equations where the linear momentum parameters are —Pi and — P2, i-e. the problem where the 
sign of all components of the linear momentum parameters has been reversed. The proof can be found in appendix A. 
It is also possible to prove that a solution where the parameters Pi, P2 = and Si and S2 has some non-zero values 
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will have the same metric as a solution with parameters —Si and —82- Therefore, configurations E and F (figure pi] ), 
will have the same total mass Mtot- The same holds for configurations B and C. 

For configurations with linear momenta, we will compare numerical results for Mtot with the predictions of special 
theory of relativity, 

Mf„™ = ^Mf + P2 + ^Mi + P^ . (39) 

Equation (^9|) ignores gravitational interaction between the holes. We also can compare Mtot of a two black hole 
configuration {M, Pi, P2) with the sum A/j*^^ of numerical Mtot of individual black holes with parameters (A/, Pi) and 

{M,P2), 

Mt*ot = Mtot {M, Pi ) + Mtot ( Af , P2 ) . (40) 

Again, Mtot ignores mutual gravitational interaction of the black holes. 

For configurations with angular momenta, we can compare numerical Mtot with 

Af,;, = MtotiM, Si) + MtotiM, S2) . (41) 

similar to the linear momentum case. 

Now we discuss configurations without linear or angular momentum and configurations with linear momentum P = 5 
and with varying separation d between the holes. For these configurations, we found that the difference between Mtot 
extrapolated to r = 00 and Mtot evaluated at r = 512 is of the order of 1%. In what follows, we present values of 
Mtot evaluated at r — 512. 



A. Two black holes with zero momenta 



Figure |lj shows Mtot as a function of d for Case A, the configuration of two black holes without linear or angular 
momentum. The exact solution for this case can be found with the semi-analytic method described in ||l5[| . As one 
can see, the total mass converges to 4.0 when the distance between the holes increases. This is expected, because 
the gravitational attraction between the holes (which also contributes to the total mass) becomes negligible as the 
distance between the holes gets larger, and the holes can be regarded as point particles. However, the gravitational 
attraction should decrease the total mass, but we see that Mtot increases or is constant with decreasing d. This 
behavior clearly shows the presence of an extra gravitational field that artificially cancels the effects of gravity on the 
total mass (in the puncture data) or adds even more energy that the amount required to cancel out the effect on the 
total mass (in the imaging data). The puncture solutions contain less junk that the imaging solutions. 



B. Two black holes with linear momenta 



Figures |l^ and ^show Mtot for cases B, C and D which contain one black hole with linear momentum P = 5 another 
with P = 0. The figures show that, in both conformal-imaging and puncture methods, Mtot is practically independent 
of orientation of P with respect to the line passing through the black hole centers. At large separations, Mtot for both 
methods does not tend to the special relativity limit Mf^^ . Instead Mtot tends to Mtot (Equation (28)). This shows 
that even at larger separations, the solutions contain a junk gravitational field. This junk field is associated with the 
junk fields of individual solution for black hole with P — 5. From the comparison at small separations we see that 
puncture solutions contain less junk than conformal imaging solutions. 

Figures ^ and 16 show Altot for cases where both black holes have linear momentum P = 5 but with different relative 
orientations. We see that Mtot for these cases depend on relative orientation of P. When Pi and P2 are antiparallel, 
Mtot is practically the same, independent of orientation of P (compare cases E, F and J). These configurations also 
have minimal masses compared to other cases at the same separations. Configuration G with parallel momenta 
gives the largest Mtot which is almost independent of d. Configurations H and I with perpendicular momenta are 
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intermediate between parallel and antiparallel cases. At large separations, Mtot seems to tend to M*^^ taking into 
account that Mtot computed at r = 512 is slightly less than Mtot at r = oo. Again, we conclude from figures ^ and 
[T6| that both conformal-imaging and puncture solutions contain junk fields. The amount of junk appears to be less in 
the puncture method. 

Figures ^ and |l8| show Mtot as a function of P for black hole configuration B with only one black hole having non-zero 
P = 5 and for configurations E and G with anti-parallel and parallel linear momenta P = 5 for the separation d = 128 
between the black holes. We see that Mtot is somewhat larger for black holes with parallel momenta than for black 
holes with anti-parallel momenta. This is the same tendency which we already saw on figures |l^ and |l^. We found 
that the values of Mtot for cases E and G are somewhat smaller that the sum of Mtot of single black hole with the 
same values of P. For case B, Mtot is very close to the sum of Mtot a Schwarzschild black hole with P = and a 
boosted black hole with linear momentum P = 5. 



C. Two black holes with angular momenta 



Figure ^ shows Mtot as a function of separation d for configurations with two black holes where one black hole has 
S = and the other one has S — 5, calculated by both conformal-imaging and puncture methods. First, Figure [l^ 
shows that Mtot does not depend on the relative orientation of S with respect to the line connecting the black hole 
throats or punctures. At large separations Mtot tends to the sum of the numerical Mtot of two single black holes with 
M = 2 and the same S — and 5* = 5, respectively. As it was for similar configurations with linear momenta (see 
figures |l3| and |lj), the amount of junk is greater in conformal-imaging solutions. 

Figure |2^ shows Mtot as a function of separation d for two black hole configurations with both black holes having 
non-zero angular momentum S = 5, calculated by both conformal-imaging and puncture methods. Again, Mtot does 
not depend on relative orientation of angular momenta. At large d, Mtot tends to the sum of Mtot for two individual 
rotating black hole with S = 5. Again, the amount of junk is greater in conformal-imaging solution. 

Figure |l] shows Mtot as a function of 5* for all two-black hole configuration with angular momenta (Figure ^l]) and 
separation d = 128, computed by both imaging and puncture methods. We found that at this separation, Mtot for all 
of these configurations practically coincide with the sum of Mtot of two individual black holes with the corresponding 
momenta. 



VI. CONCLUSIONS 



We calculated conformal initial value data for single and binary black hole configurations on a high-resolution, adap- 
tively refined Cartesian mesh using both conformal-imaging and conformal-puncture methods. Adaptive mesh re- 
finement approach allowed us to obtain accurate three-dimensional numerical solutions both near and far away from 
black hole for a wide range of black hole configurations with various linear and angular momenta and various black 
hole separations. The equivalent uniform-grid resolution obtained in the calculations was 65536'^. 

Determination of the total mass Mtot of configurations required an integration over a surface at infinity (|l2|). Using 
numerical values of Mtot obtained by integration over a sphere at large finite radii r < SOOrg, where rg — 2M/c^ 
allowed us to obtain M^ot (oo) by extrapolation. For configurations with large linear momenta, the difference between 
Mtot(oo) and Mtot(500rg) was as large as ~ Q.lMtot- The difference was substantially smaller for configurations with 
large angular momenta. Comparison with known analytic solutions was then performed using asymptotic values of 
Mtotioo). 

For a single rotating black hole, the mass of the calculated apparent horizon Al^h is less than the irreducible mass of 
a Kerr black hole with the same Mtot- For a single boosted black hole, the calculated Mah is less than that of the 
corresponding boosted Schwarzschild black hole. This shows that the junk field is present in the conformal solutions. 
Comparison of the spatial distribution of local invariants for both rotating and boosted black hole shows that junk 
field is present at large distances from the holes. For a rotating black hole, the shape of the numerical apparent 
horizon is much more spherical than that of a Kerr black hole. It appears the conformal approach generates more 
spherically-symmetric solutions. 
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For two-black hole configurations, with increasing separation d between the black holes the value of Mtot{d) tends to 
the sum of Mtot of two individual single-black holes numerical solutions containing their corresponding junk fields. If 
only one of the black holes has a non-zero linear or angular momentum, the value of Mtot is practically independent of 
the orientation of the momentum. For black holed both having non-zero angular momentum, Mtot is also practically 
independent of the orientation of the momenta. For black hole both having non-zero linear momenta, the value of 
Mtot is maximal for parallel momenta and minimal for anti-parallel momenta. In general, the conformal-puncture 
method appears to generate less junk gravitational field than the conformal-imaging method. 

This work confirms that conformal initial conditions are not suitable as initial conditions for astrophysical black hole 
collisions since they contain junk gravitational field. They are useful as initial conditions for testing time- integration 
schemes for black hole collisions because black holes can be places at small initial separation. Our adaptive mesh 
refinement approach allows to construct these initial data on a Cartesian mesh with a high accuracy. We plane to use 
similar adaptive mesh refinement for the integration of these initial data in time. 
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APPENDIX A: APPENDIX: PROOF 



Proposition: Let a configuration of two black holes be given. Suppose that the linear momentum parameters of 
the problem are PI and P2, and that the angular momentum parameters for both holes are the null-vector. Then 
the metric that is a solution to the constraint equations with these parameters is the same as the metric of that is 
the solution to the constraint equations where the linear momentum parameters are —PI and — P2, i.e. the problem 
where the sign of all components of the linear momentum parameters has been reversed. 

Proof: The energy constraint is given by: 

V^^+i^-^AyXy =0 , (Al) 



As we can see, the extrinsic curvature enters only as KijK^K Thus, to prove the proposition, we need to prove that 
this term does not depend on the sign of the linear momentum vectors. In the puncture formulation the extrinsic 
curvature tensor of the configuration PI, P2 is given by: 

K^i = — (PVnP + PVnV - (f^ - nVnP) Pl^nlA + 
2rl ^ \ / / 

^ (P2*n2^ + P2%2* - {f^ - n2'n2^) P2^n2i) (A2) 

Trivial manipulations give that: 
q / 1 

KijK'^ = 4 ( "'^^ + "-^^ ~ ^ "■^'^ ^'^'^ 

+nli nlc nV nP Pl^ Pl'^ - 2 nh nlc PI' Pl'^ + 3 {nlcf {Pl^'f- 
2 nli {nlcf nV {Pl^f + nlj nV PP (Plj + nU nl^ PI"") + 
nljTiV {Pli + nlinl^Pl") {PV- + nl^nr PV)) 
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+ ^^2^2^ (nl^' n2j PV P2, + nl' n2j Pl^ P2, + nl^ nV n2j Pl^ P2, 
~2 nlc n2' PI"" P2i + nV n2i PV P2j + nl^ n2, PV P2j + 
nlc nV nV n2, P2j + nl, n2^ Plj P2' - 2 nl^ n2, PI'' P2'+ 
nlj n2^ {Pli + nl, nlc PI"") P^' + nli n2' Plj P2^ + 

nlj n2' {PI, + nli nlc PI") P'2^ - 2 nV n2c PI, P2'= + nl, n2c n2' n2^ Plj P2'=- 

2 nU n2c PV P2^ + nP n2i n2j n2c Pr P2^ + nV n2i n2j n2c PP' P2'= 

+6 nlc n2c PV P2'= - 2 nl, nlc nV n2c PI" P2^ + nlc nV n2, n2j n2c PI" P2'=- 

2 nlc n2, n2c n2' PI' P2^ + nlj n2c n2' nV [PI, + nl, nlc PI") P2") 

+ ^ {n2j n2^ P2, P2* + n2i n2^ P2j P2' + n2j n2' P2, P2' + n2, n2' P2j P2^ 

-2 n2c nT P2i P2'= + n2j n2c n2' n2^ P2, P2'= + n2, n2c nT nV P2j P2^ 
-2 n2i n2c P2* P2^ + n2i n2j n2c nV PV P2^ + n2, n2j n2c n2' PV P2=+ 

3 (n2c)^ (P2'=)2 - 2 n2, {n2c)^ nT [P2''f + n2i n2j (n2c)^ n2' nV {P2''f)) 

What's interesting about this equation is that each term have the product of two components of the momentum 
parameter vectors. Thus, if we change the sign on all components of the momentum parameter vectors, it will not 
change anything in the value of the KijK^^ term, and thus the metric is conserved under such a sign change, at least 
in the puncture formulation. In the Cook formulation, there are more terms in the expression for Kij and K^^ but 
each of these terms is a product of some factor that depends on the position where we wish to find Kij and of the 
extrinsic curvature tensor, evaluated at some different position. This means, that we multiply two terms where the 
normal-vector and the the radius that enters in the expression for the extrinsic curvature tensor are different, but the 
parameters PI and P2 does not change. Thus each term in the product KijK^^ will contain a product of two terms of 
the momentum vector parameters, and thus the whole expression is conserved under the sign change of Pi, P2. QED. 
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FIG. 1. Asymptotic total mass, Mtot, and total angular momentum, Stot, for case M=2 S=40. The diamonds are the 
numerical values of Mtot, Stot found from integrating over a sphere at radius r. The line is a least square fit of the data to a 
second order polynomial. 




Mass vs. S, 1 black hole, imaging data 
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FIG. 3. Total mass, A4tot, apparent horizon mass, Mah, and irreducible mass, Mirr, as a function of S for a rotating black 
hole. Mtot and Mah are calculated numerically, and Mirr is the irreducible mass of a Kerr black hole with the same overall 
mass and angular momentum as the numerical black hole. 
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FIG. 5. 7i (eq. |T|), as a function of isotropic coordinate y (which coincides with the equatorial axis) for rotating numerical 
and Kerr black holes. 
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FIG. 6. Relative comparison of h (eq. |l]) at large distances for rotating holes. 
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FIG. 7. Asymptotic total mass, Mtot, and total linear momentum, Ptot, for case M=2 P=40. The diamonds are the numerical 
values of Mtot, Ptot found from integrating over a sphere at radius r. The line is a least square fit of the data to a second order 
polynomial. 
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FIG. 8. Total mass, Mtot, apparent horizon mass, Mah, and boo sted mass, Mioosted, as a function of P for a boosted black 
hole. Mtot and Mah are calculated numerically, and Mtoosted = \J M^^^ — jg the horizon mass of a boosted Schwarzschild 
black hole with the same overall mass and linear momentum as the numerical black hole. 
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FIG. 9. Ii (eq. p^ ), as a function of isotropic coordinate x (which coincides with the equatorial axis) for moving numerical 
and boosted Schwarzschild black holes. 
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Mtot vs. d, Case A, imaging and puncture 
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FIG. 12. Total mass, Mtot as a function of separation, d, for two black holes with configuration A (figure pT| ) 
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FIG. 13. Total mass, Mtot, as a lunction ol separation, d, for two black holes with configurations B, C, D (figure |ll|) with 
linear momentum, imaging method 
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FIG. 14. Total mass, A4tot, as a function of separation, d, for two black holes with configurations B, C, D (figure |l^ with 
linear momentum, puncture method 
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Mtot vs. d, Case E+F+G+H+I+J, imaging 
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Mtot vs. d, Case E+F+G+H+I+J, puncture 
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Mtot vs. P, Case B+E+G+lbh, imaging 
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Mtot vs. P, Case B+E+G+lbh, puncture 
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Mtot vs. d, Case B + C+D, imaging + puncture 
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FIG. 19. Total mass, Mtot, as a function of separation, d, for two black holes with configurations B, C, D (figure |llj) with 
angular momentum 
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Mtot vs. d, Case E+F+G+H+I+J, imaging + puncture 



M 



tot 



5.3 



5.2 



5.1 




CaseE+F,img 
Case J,inig 
Case G,img 
Case H+I, img 
Case E+F,punc 
Case J,punc 
Case G,punc 
Case H+I, punc 
Mtot* img 
Mtot* punc 



20 



40 



60 



80 



100 



120 



27 



FIG. 20. Total mass, Mtot, as a function of separation, d, for two black holes with configurations E, F, G, H, I, J (figure 
with angular momentum 
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FIG. 21. Total mass, Mtot as a function of angular momentum, 5", for separation d = 128 for two black holes with configu- 
rations B, E, G (figure 0) 
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